Theoretical investigation of magnetoelectric effects in Ba2CoGe207 
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A joint theoretical approach, combining macroscopic symmetry analysis with microscopic methods 
(density functional theory and model cluster Hamiltonian) , is employed to shed light on magneto- 
electricity in Ba2CoGe2 07. We show that the recently reported experimental trend of polarization 
guided by magnetic field (H. Murakami et ah, Phys. Rev. Lett. 105, 137202 (2010)) can be pre- 
dicted on the basis of phenomenological Landau theory. From the microscopic side, Ba2CoGe2 07 
emerges as a prototype of a class of magnet oelectrics, where the cross coupling between magnetic 
and dipolar degrees of freedom needs, as main ingredients, the on-site spin-orbit coupling and the 
spin-dependent Op- Co d hybridization, along with structural constraints related to the non- 
centrosymmetric structural symmetry and the peculiar configuration of C0O4 tetrahedrons. 

PACS numbers: Vahd PACS appear here 



The magnetoelectric (ME) coupling between the mag- 
netic order parameter (M) and the ferroelectric one (P) 
has been well studied in recent years in the context of 
multiferroic oxides, where two or more primary ferroic 
phases coexist in the same system. [l| Most of the re- 
search focused onto the linear ME effect, as expressed by 
the aijHiEj cross-coupling term, and observed, e.g., in 
the prototypical Cr203.[l] Microscopically, the ME effect 
can be ascribed to i) the inverse Dzyaloshinskii-Moriya 
mechanism Q (or the equivalent spin-current mechanism 
[1]) showing P (X Cij x {Si x Sj) between neighbor- 
ing spins connected by a vector e^j, and ii) the inverse 
Goodenough-Kanamori (or exchange-striction) mecha- 
nism [El showing P (X Jij{Si-Sj) with exchange inte- 
gral Jij. A third mechanism has been recently proposed, 
namely iii) the spin-dependent p-d hybridization [6[ 0, 
where the spin-orbit-coupling (SOC) "asymmetrizes" the 
p-d hybridization between the transition metal (TM) and 
the surrounding ligands, inducing an electric polariza- 
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tion P oc ' ^jJ ^j^ where labels the vec- 

tors connecting the TM to the ligand ions. Since the 
third mechanism can concomitantly occur with the first 
or second one (in non-collinear spin structures), it is in 
general difficult to identify each of these mechanisms. 
However, it has been recently reported that the third 
mechanism can be responsible for the polarization ob- 
served in Ba2CoGe207(BCGO), where two neighboring 
Co spins are aligned in an antiferromagnetic (AFM) con- 
figuration. Incidentally, we have discussed a related ME 
mechanism in Fe304, where the crystal structure with 
Cc space group doesn't have inversion symmetry, but 
the c-gHde symmetry prohibits polarization along the y 
direction, [s^ Taking into account SOC and ferrimagnetic 
spin order, i.e. considering the magnetic group, the sym- 
metry is broken and a small spin-dependent contribution 
to polarization arises. In this letter, we apply a similar 
theoretical analysis to BCGO: we first perform a sym- 
metry analysis, then show the DFT results on the ME 
effect, and further confirm the microscopic mechanims in 
a model Hamiltonian approach. 



FIG. 1: Crystal structure in the ab plane; Co and Ge ions 
(both located in O4 tetrahedron) lie in c=0 planes whereas 
Ba ion lies in c=l/2 planes. Co spin configurations: (i), (ii) 
collinear AFM and (iii) non-collinear spin-canted under ap- 
plied H//no. 



Ba2CoGe207 (melilite) shows tetragonal non- 
centrosymmetric (but non-polar) P42im (#113) 
structure with two Co sites. Col at (0,0,0) and Co2 at 
(1/2,1/2,1/2) sites, as shown in Fig. E Below TAr=6.7K, 
the magnetic structure shows collinear AFM spins lying 
in the ab plane. Experimentally, the electric polarization 
along the c axis, Pc, is measured even when ij"=0.[9[ fiol 
Additionally, P develops finite components along any 
direction when H is applied, modulated by the direction 
and the size of the magnetic field. 

Symmetry Analysis. — In order to characterize the pe- 
culiar ME effect, we first briefly discuss the group theory 
analysis and its implications in the framework of Landau 
theory of phase transitions. [11] In the parent PA2iml' 
magnetic space group with eight symmetry operations 
C2{z)-> 2'S'4, 2C2{x,y)-> 2cr(i} plus time-reversal (1'), the 
magnetic order leads to a lowered symmetry. We define 
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TABLE I: Matrices of the generators of space group 
P42iml' in the representations spanned by F, A and P. 
The group elements denote the identity, vr-rotation, 7r/2- 
rotoinversion, screw C2y+(^^0) and time- reversal. Labels of 
irreducible representation (IR) are taken from the ISODIS- 
TORT program. [12 
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the order parameters F = S'i+S'2 and A = S'l — S'2 as the 
ferroraagnetic (FM) and AFM combination of Col and 
Co2 spins, respectively. Using the transformation rules 
given in Table HI we express the thermodynamic free en- 
ergies in terms of all the possible ME coupling terms of 
the form P ■ M'^ which are invariant under symmetry 
operations: 

i^ME = CAPcAaAb + C^PcFaFi, + CA^PMaFa - A^Ffo) 
-^Ci{PaAaFc - PbAbFc) + C2{PaAcFa - PhAcF^), (1) 

while the dielectric energy is Fde = -PV2x, where 
ca^cf^ci^C2^caf and x (henceforth set as 1) are con- 
stants. P is then evaluated at the minima of F = 
Fme + i^DE, reading 

Pa = CiAaFc + C2AcFa, Pb = -CiA^Fc - C2AcFb, 

Pc = CAAaAb + CpFaFb + CAF{AaFa - A5F5). (2) 

Note that only the first term of Pc originates purely from 
the AFM order, explaining the observed spontaneous P, 
whereas other components are allowed only in the pres- 
ence of the FM order parameter. 

Hereafter we focus on the Pc behavior assuming a 
canted AFM configuration under an applied magnetic 
field, i.e., we first simultaneously counter-clock- wise ro- 
tate two antiparallel Co spins in the ab plane with the 
angle (f) from the a axis, then we cant spins by an an- 
gle as depicted in Fig. [H Accordingly, we set Si = 
S ( cos{(j) + sin( + 0) and 82= S{ — cos{(j) — 
(j)')^ — sin((/) — 0), ending up with 

Pc{^) = sin 2(/)(ca cos^ cj)' — cp sin^ cj)' — cap sin 2(j)') 
= 2aS'^ sin 2(/) cos(2(/)' - f3) , (3) 

where = c^p + (ca + cf)^/2 and tan/3 = — (ca + 
cf)/2caf. By neglecting the canting angle Eq. ([3]) 
perfectly reproduces the experimentally observed depen- 
dence of polarization on the spin angle, being Pc ex sin 2(j) 
at T=2K and iJ=lT.[l^ A spontaneous Pc can be there- 
fore induced in the Ano (Ai_io) order but not in the 
Aioo (Aoio) one (cfr. Fig. [Tj). Analogously to the case of 
magnetite 0, the non-magnetic group lacks the inversion 



symmetry, but the symmetries which prohibit Pc (e.g. 
C2y rotation) are broken by the Ano magnetic order. 
Starting from the Ano order, further symmetry reduc- 
tion occurs by applying an external H. Indeed, Ano 
order shows 2'^ point group, which allows non-zero ais, 
<^23, <^3i, <^32 linear ME components [l^ in such a way 
that Pa and P5 can be induced by applying H^. Fi- 
nally, Eq. (|3|) at fixed (j) gives the simple ^'-dependence 
Pc{(l)') oc acos{2(j)^ — where the phase shift depends 
on the non-zero cap coefficient. 

DFT analysis — In order to quantitatively confirm the 
ME behavior and to investigate its microscopic mecha- 
nism, we performed DFT calculations using 
with GGA-PBE potential (we checked our results by us- 
ing also GGA-\-U\TE\ potential with U=3 or 5 eV for 
Co-d state). Due to the lack of experimental information 
on structural parameters, we considered the Ca2CoSi207 
structuredll and optimize it by substituting atoms (Ca 
^ Ba, Si o> Ge) without SOC. The optimized struc- 
tures shows a=6=8.28A and c=5.58A, and the tilting 
angle of C0O4 tetrahedron given by Ai:=23.9°, consistent 
with experimental values, a=6=8.41A and c=5.54A0 
and K.=2r.^ 

TABLE II: Magnetic anisotropy energy (MAE) (meV/Co) 
obtained by comparing the total energy with different spin 
directions under SOC and for different values of U in the 
GGA+t/ scheme. Spin and orbital moment (/xb) are also 
reported for *S'//(100). In the rightmost column we report 
the calculated Pc {iiC/m^) for S (L)//110 with fixed atomic 
structure. 
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In the C0O4 tetrahedra, the Co^+ ion shows orbital- 
quenched e'^^t2^ge'^^t2g states, which causes a very small 
magnetic anisotropy, as shown in Table |TT1 The ob- 
served magnetically easy ab plane and hard c axis are 
consistent with experimental report (S//010 from Neu- 
tron diffraction! 17]). The small MAE in the ab plane 
explains why the spins easily follow an applied H: even 
under a small magnetic field, the spins flop to be per- 
pendicular to H and then cant in order to reduce the 
Zeeman energy. 

Imposing the collinear AFM configuration, we simul- 
taneously rotate the Co spins in the ab plane. We eval- 
uated the ME effect as the change of P (calculated by 
Berry phase method [HI) induced by the rotation of M 
with respect to the crystalline axes, in the fixed non-polar 
crystal structure. In Fig. [2] (a) we show Pc as a func- 
tion of the spin-rotation angle 0, consistent with both 
experiments and the previously discussed Landau analy- 
sis, being Pc oc sin 20. The calculated polarization, which 
displays a maximum value Pc = 12.1 jiC /w? at Ai_io, 
originates here from a purely electronic contribution via 
SOC, and is further enhanced when atomic internal co- 
ordinates are optimized in the Ai_io configuration, as 
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FIG. 2: (a) DFT results for Pc as a function of the collinear 
spin angle (j) in the ab plane, fitted to /((/>) = —a sin 20, with 
a=12.7 (solid line). Spin configurations in the ab plane are 
shown by arrow, (b) DFT results for Pc as a function of the 
noncollinear spin-canting angle 0^ in the ab plane, fitted to 
= a cos 2(0' - 6) + c, with a=17.9, 6=22.7 and c=0.06 
(solid line). Inset: Pc as a function of magnetic field. 

discussed later. 

TABLE III: Pc at different canting angle 0', calculated in 
the fixed non-polar structure (first line) and with optimized 
internal atomic coordinates (second line) . Experimentally, 
the maximum of Pc is about 120/iC/m^ [10]. 



Pc (/iC/m") 


=0° 30° 90° 


fixed structure 
opt. structure 


12.7 17.6 -12.7 

39.8 57.7 -38.9 



We look then at the spin-canting effect induced by 
an applied field Hhq. In Fig. [2] (b) we show the 
change in Pc induced by artificially canting the spins 
by an angle starting from the Ai_io AFM configu- 
ration. In agreement with the Landau theory analysis, 
Pc evolves as cos 2(0' — 22.7°) + const., displaying a peak 

(j)' ^ Hi. We also evaluate the evolution of Pc as a func- 
tion of the applied H (cfr. inset in Fig. [2](b)), assum- 
ing the experimentally measured magnetic susceptibility 
X = M/Hiio -0.25 /ie/T per Co[10]. Although the 
trend of P shows good agreement with experiments [lo|, 
its size is one order of magnitude smaller. This devia- 
tion is reduced when the atomic structure is optimized 
in the canted-AFM configuration, as shown in Table 
mil This means that ferroelectricity is strongly coupled, 



through magnetism, with lattice distortions in a sort of 
magnetically-induced piezoelectric effect. The PH curve 
— first increasing with H and then decreasing and chang- 
ing its sign — denotes an atypical ME trend. Although 
one may assume a non-linear ME coupling coming from 
high-order terms in the free energy, the nontrivial evo- 
lution of P actually arises from the composition of the 
three P ■ M'^ terms appearing in Eq. ([1]). 

Single-site SOC induced ME effect — The microscopic 
origin of P in BCGO can be easily explained in terms 
of a cluster Hamiltonian for a single C0O4 tetrahe- 
dron, which allows to further clarify the role of the lo- 
cal SOC in the spin-dependent p-d hybridization mech- 
anism. Neglecting contributions from the energetically 
deeper majority-spin states, the Hamiltonian consists 
of four terms, H = Hd -\- Hp -\- Hpd + i^soc, where 
Hd = ^T^a^Ua and Hp = £p Ez,/5^I/3^/,/3 account 
for the local energies on Co and O sites (with Sp = 
the energy reference and A = Sd — £p), which hy- 
bridize through Hpd = ^o^l^M^ocPi,^ + ^•^)- ^^^^ 
a and {3 refer to the d = xy, yz^ zx^ x^-y'^^ Sz^-r^ and 
p = x^y^z orbitals involved, whereas / = 1,...,4 la- 
bels the four oxygens surrounding the Co ion, located 
at Ri = (1,1,-1), (-1,-1,-1), (1,-1,1) and (-1,1,1) 
in the local reference system with Co in the origin. 
The hybridization matrix Va/3i depends on the d and 
p orbitals involved (with a or tt bonding) and on the 
relative positions of the ions; we adopted the Slater- 
Koster parametrization[19j , assuming tpdg- = 1.3 eV, 
A = 5.5el/[20| and tpd^r = -0.45 tpda^- The last 

term is i^soc = A " ^W)^ida'^ where the 

matrix elements can be expressed as a function of the 
polar and azimuthal angles (^, 0) defining a local ref- 
erence for the spin-quantization axis[22l]. We assume 
A = 0.021eV, the free Co ion spin-orbit coupling value. 
SOC-induced mixing of the local d levels lifts the de- 
generacies in the Cg and t2g manifolds and implies dif- 
ferent hybridizations with the ligand oxygens, that may 
ultimately induce a local dipole moment. We evaluated 
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-90 90 180 

spin angle in ab plane (|) 

FIG. 3: Model results, (bottom) Electron density of upper- 
and lower- lying oxygens, po, in the tetrahedron model as a 
function of the azimuthal angle (f) at = 90°. (top) Orbital 
occupancy on Co, p^Q, as a function of 0. 
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then the local occupancies as po(0 = (Pi ^Pi (b) 
PCo — i^a^a) ^ fuuctiou of the aziuiuthal angle 
0, i.e. rotating the spin in the ab plane. As shown in Fig. 
[31 we can distinguish between lower-lying (Oi, O2) and 
upper-lying (O3, O4) oxygens, with pup,io ^ ±sin(2(/)). 
Then a local dipole p = (e/4) J]]^ po(0 Ri de- 
velop only along c, proportional to the charge difference 
Apo = pup—pioi i-e. Pc 2sin2(/), in excellent agreement 
with the predicted functional form P ex "^^^jiSi • e^)^ e'j 
0, S [IqI- Furthermore, we can estimate the (i— orbital 
mixing on the Co site by looking at the orbital occupan- 
cies, shown in Fig. O Even if the two occupied states 
have prevalent dx2_y2^ds^2_j.2 characters, a small mix- 
ing occurs via SOC with (mostly) dyz^dzx orbitals, be- 
ing pyz (X cos^ (/), pzx ^ sin^ ^, i.e. the most occupied 
is the one perpendicular to the spin-quantization axis. 




non-bonding asymmetric-bonding 



FIG. 4: Bonding nature of (i-orbital in O4 tetrahedron and 
induced local polarization. Via SOC, asymmetrically bonded 
orbital states are mixed with non-bonding occupied states. 
Inset: DFT-calculated energy levels of orbital states. Possible 
SOC mixing in minority spin states are shown with energy 
difference A^. 

TABLE IV: DFT-calculated 3d orbital-decomposed occu- 
pancy (in percentage, with spin states summed up) with 
different SOC enhancement factors A (0= without SOC, 
xl=with standard SOC, and xlO=with the SOC term 10 
times artificially enhanced) for different S directions in local 
xyz frame. 



A M 


'2 '2 '2 '2 

xy yz zx Sz -r x -y 





50.0 50.0 50.0 100.0 100.0 


xl S//x 
xl S//y 


50.0 50.2 50.0 99.9 99.9 
50.0 50.0 50.2 99.9 99.9 


XlO S//x 
XlO S//y 


49.5 61.7 49.7 92.8 96.3 
49.5 49.7 61.7 92.8 96.3 



These findings are in excellent agreement with DFT cal- 
culations, as shown in Fig. |4] and in Table |lVl where 
the hierarchy of d-orbital occupancies at selected values 
of the spin direction is confirmed. Such a mixing of lo- 
cal d-levels nicely explains why Pc size is maximum at 
(j) = ±45°, when the Co spin is parallel either to the 
upper- or to the lower-lying oxygen bond; indeed, as pic- 
torially shown in Fig. 21 the composition of yz and zx 



orbitals has an asymmetric bonding nature in the tetra- 
hedron, i. e. non-bonding with upper ligands and bond- 
ing with lower ligands or vice versa. 

We can further estimate from our model the Pc de- 
pendence on the spin angle in the system by considering 
two C0O4 tetrahedra tilted by k,. In the AFM collinear 
configuration we find Pc oc pd k) -\- Pc2 (^ + tt — /^) = 
2 cos ti sin 2(j). Analogously, in order to mimick the effect 
of the external i^no, we can define the canting angles 
diS = (/) + — 7r/4 and (j)'2 = —(j) + + 37r/4, find- 
ing Pc{4>') oc cos2{(j)' — n)^ in excellent agreement with 
experiments and DFT results. 

Conclusions — 

We shed light on the mechanism underlying pecu- 
liar magnetoelectric effects in Ba2CoGe207, by combin- 
ing different theoretical approaches and explicitely tak- 
ing into account the microscopic atomic arrangement 
and symmetries of the compound. Our Landau phe- 
nomenological theory shows that: i) On top of non- 
centrosymmetric non-polar P42im symmetry in the non- 
magnetic crystal structure, a collinear antiferromagnetic 
spin configuration with in-plane spins allows an elec- 
tric polarization along the z axis, ii) Upon applying an 
external magnetic field, the induced non-collinear spin- 
canting well reproduces the experimentally observed pe- 
cuHar trend of polarization related to the tilting angle 
between C0O4 tetrahedrons. In order to have quantita- 
tive estimates, we perform relativistic ab-initio calcula- 
tions and highlight the delicate interplay between orbital 
occupation and local magnetic anisotropy, resulting in 
an excellent match with available experiments. Further- 
more, as a proof that the microscopic origin of magneto- 
electricity is based on two relevant ingredients (i.e. the 
anisotropic p-d hybridization between Co and O states 
and the on-site spin-orbit coupling), we built a realistic 
tight-binding model, taking into account C0O4 tetrahe- 
dron and the crystal field splitting, that is sufficient to 
nicely explain magnetoelectric effects and put forward 
Ba2CoGe207 as a prototype of the class of materials 
where the interplay between magnetism and ferroelec- 
tricity is based on spin-dependent p — d hybridization, 
as recently suggested in the literatures based on two-ion 
cluster model. ^0|, HI 

During completion of the work, we became aware of 
a similar symmetry analysis [2^ performed for BCGO 
by Toledano et al. However, their focus is on toroidal 
moments, whereas ours is on the combination between 
single-ion anisotropy and p — d hybridization (derived 
from density functional and tight-binding models) as 
main microscopic mechanism driving magnetoelectricity. 
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